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Abstract 

We consider global schemes with 7? 2 -estimates for the multivariate 
Burgers equation and the incompressible Navier-Stokes equation in its 
Leray projection form. We extend the definition of a global scheme 
and the global L 2 -theory of the multivariate Burgers equation to the 
incompressible Navier-Stokes equation, where we consider estimates 
of the gradient of the pressure in its Leray projection form in H 1 at 
each each time step of the scheme. We use the fact that the data in 
our scheme are in H 2 at each time step such that the data for the 
Poisson equation related to the elimination of pressure are in H 1 at 
each approximation step. The growth of the velocity components is 
linear with respect to the if 2 -norm on a time scale of time step size 
pi ~ j at each time step I > 1. The time-step size is 'minimal' in order 
to make the scheme global while at the same time it provides us with 
uniform bounds for the first order coefficients of linear approximations 
of the local solutions at each time step. 
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1 Introduction 



In [8] we considered a constructive global scheme for the multivariate Burgers 
equation 
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u(0,.) 



on a domain [0, oo) 



l n , and where the viscosity constant v is strictly pos- 
itive, i.e., v > 0. We assumed hi € H s = H s (W 1 ) for s € M and do this 
henceforth in this paper. Furthermore, in order to have a uniform represen- 
tation of the following argument we assume that n > 3. The differences to 



the Cauchy problem for the incompressible Navier-Stokes equation, i.e., 
* % - i/Av + (v ■ V)v = -Vp, 



V-v = 0, t>0, xeM. n , 
{ v(0,.)=h, 



(2) 



are the source term in form of the negative gradient of the pressure, and the 
incompressibility condition. A classical solution with velocity components 
Vi, 1 < i < n, and such that Vi(t, .) € C 2 n H 2 for all t solves 



dt v luj=l -Qxf + L>j=l u dxj 

{ v(0,.)=h, 

for 1 < i < n, and vice versa. In this case the Poisson equation 

-Ap = E?, fe =iGfe^. 
has the well-defined solution 
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is the Poisson kernel. Here, we mention the case n = 2 and invite the reader 
to apply the following argument in the case n = 2. However, since this case 
is different from all cases n > 3 we assumed n > 3 above and stay with this 
assumption. We also mention that |.| denotes the Euclidean norm and u> n 
denotes the area of the unit re-sphere. In equation ([3]) we use the fact that 
the formula @ has a well defined gradient such that 



Vp(t,x) = / VK n (x-y) V 

jRn j,k=l 
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dxj dxk 
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This is the term we need to control in order to extend the scheme for the 
multivariate Burgers equation in [8] to the present situation. Note that for 
n = 3 the functions 
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are not integrable outside a ball B^(0) := {x||x| < R} with origin and 
radius R > 0, i.e., are not in L 1 (W 1 \ B^(0)). It is useful to note that we 
have local integrability in the Li-sense of the first partial derivatives of the 
kernel K. For n > 3 we also have square integrability since 
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Integrability in the L 1 -sense of functionals as in ([7]) is harder to achieve, even 
if we invoke partial derivatives using some information about the gradient 
of velocity. Note that even if we consider the second derivative we of the 
kernel, then we get 
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if i = j, 
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and these functions are surely square integrable outside a ball around the 
origin for n > 3 but not intergable in the L 1 sense outside a ball around 
zero. Now the gradient of the pressure in ([5]) contains a source of a Poisson 
equation in the form 

^ f dvk dvj 
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dxj dxk 
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If our approximations v?'" 1 ' 1 of the velocity components in (jlip at each 
time step / are in H 2 and have a uniform bound for all iteration step numbers 
m of the approximation such that we have a bound C + IC that is linear with 
respect to the time step number I, then we have the possibility to extend 
L 2 -theory for constructive scheme of the multivariate Burgers equation [8] 
to the incompressible Navier-Stokes equation. In this case we have for all 
approximations the pointwise relation 
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and the linear bound leads to the estimates 
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For the estimation of convolutions we use the generalized Young inequality, 
i.e., the fact that for 1 < p, q, r < oo 

f £ L q and g (z L p -> / * g G L r , if - + - = 1 + -, (14) 

p q r 

and 

\f*9W < \f\Lp\g\Li (15) 
for a convolution f * g. So (|13p together with (j!4[) and (|15p for r = 2 and 
p = 1 (estimate of the source), and q = 2 (L 2 -estimate of the source outside a 
ball around zero) lead us close to an L 2 -bound for the approximations of the 
gradient of the pressure. In order to show that our scheme is global next to a 
specific time scale we use properties of the convolutions (shifting derivatives) 
and the fact that the velocity approximations are in if 2 . In order to achieve 
estimates in different function spaces we can extend this idea observing 
that the source in (|lip contains derivatives which may be shifted to the 
kernel (by some elementary estimates and using partial integration) outside 
a ball of origin zero while inside this ball we may use the local regularity 
of the function in ([5])- or its counterparts in the scheme. Now in a scheme 
we may consider local linearizations in which an approximation function 
/>/>>^ 1 < j < n follows (with respect to an iteration) an approximation 



v 



function v 9 -' p ' 1 , 1 < j < n on a domain DJ := [I — 1,1) x M. n with step size 
pi ~ j and transformed time r = p\T (similar as in [8]). We then have to 
find a linear bound on the growth of the velocity function components, i.e., a 
bound for differences of form (using Einstein notation for usual derivatives) 

v{' p '\t,x) -vf' p ' (t,x) = 

= -Pi kn Y?j=i ifj ~ 9j) 0, 2/)%— 0, y)r l f (r, x; s, y)dyds+ 
Si- 1 Pi Jr" /r~ K n,i (z ~ y ) X (16) 
E",fe=l (fk,j+9k,j) (s,y) 



((fj,k ~ 9j,k) (s,y)) jT t f (r,x;s,z)dydzds, 



and where TK is the fundamental solution of the linear (scalar!) parabolic 



equation 



For n = 3 the natural bound is with respect to a H 2 -norm, and we shall 
see that can get a bound is with respect to the i^ m -norm for m arbitrary as 
well. A linear bound is sufficient by the choice of step size pi ~ \. 
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This is the program of this paper. In the next section we reconsider the 
global scheme for the multivariate Burgers equation. Then in the last section 
we extend our considerations to the incompressible Navier-Stokes equation. 



2 The global scheme for the multivariate Burgers 
equation 

Let us first reconsider the constructive approach to global existence for the 
Cauchy problem of the multidimensional Burgers equation in ([jrj. The fol- 
lowing considerations apply also to a class of initial- value boundary problems 
on different domains [0, oo) x Q with Q C R n . We assume smooth initial data 
h = (h\, ■ ■ ■ ,h n ) T where for all 1 < i < n the functions hi have polynomial 
decay at infinity- this is essentially equivalent to the assumption hi € H s 
for all 1 < i < n and s S K. We review some ideas in [8] from a slightly 
different point of view. Only ideas of proofs are given. We provide the 
details when we consider the extension to the incompressible Navier-Stokes 
equation. The idea of the constructive approach to global existence and 
regularity for this Cauchy problem and related Cauchy problems is to set 
up an iteration scheme with respect to time r = pit at each time step I > 1 
and with a time step size 

Pi ~ y, (18) 
and such that at each time step an equation 



du P ' . . \-^n d 2 u p / p.l du p / 



3=1 -g£- ~ PI £j=l 



U ■ 



•3 dx 3 



(19) 



u P' l (l-l,.) =uP> l - l (l-l,.) 



is solved on the domain [I — 1,1] x IR n . The solution is constructed via a 
functional series 



^ + ^5<' fc+1 '',l<z<n, (20) 



where uf' 1 '' solves 



du P ' 1,1 ( r-vn d 2 u P,1 ' L ^ n p,l—ln i \du P ' 1,1 \ 1 . . , 

-B7- = Pl l I/ Ei= 1 -^--Ej=i«i ( l - l i-)^fc-y ±<i<n, 
uP^ l (l - 1,.) = u'- 1 ^ - 1,.), 



(21) 
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and Su^' k ' 1 , 1 < i < n solves 

-piE,«^), (22) 
^ «y u p - fc+1 ' , (i-i,.) = 0, 

and where Su^ ,k+1 ' = u? ,k+1 ' 1 — uf' fc,/ for k > 1, and Ji^' 1 ' 1 = it?' 1 ' 1 — /ij 
at the first time step. We are interested in Cauchy problems of physical 
interest, and for this purposes it suffices to assume that 

hi£H s = H s (R n ) (23) 

for all 1 < i < n and s € M. Although the representations of the classical 
solutions of the members of the functional series (<5ti?' fc '') in terms of funda- 
mental solutions are not convolutions we may estmate them by convolutions 
using Gaussian a priori estimates of fundamental solutions. In order to do 
this we may use the Levy expansion of fundamental solutions as we did in [5] 
or we may use properties of the adjoint. For the estimation of convolutions 
we use the generalized Young inequality (|14p and (|15p above. This leads to 

Lemma 2.1. Let u? ,l ~ 1 (l — 1, .) 6 H s for all 1 < i < n and for any given 
s£l. Then 

uf'\r, .) G H s and 5uf\r, .) G H s (24) 
uniformly in r and for all 1 < i < n and all s£K. Furthermore, 

oo 

uf(T,.):=uf> l (T,.)+J2K k ' l £HS. (25) 

fc=i 

Proof. We consider the case n = 3 (the generalisation to n > 3 is similar). 
This is proved for H 2 first with classical representations of the fundamental 
solution. The latter exists for n = 3 since first order coefficients for the 
equation for u^ 1 ' 1 are of form u p ^~ X G H 2 and hence Holder continuous. 
Inductively with the substep number k > 1 the same is true for the first order 
coefficients of the equation for 8u^ k+1,1 . Then we may represent the solutions 
for U?' 1 ' 1 and for 5u? ,k+1,1 in terms of the fundamental solutions. Using the 
adjoint and partial integration we may represent the second derivatives of 
the solution by a representation which contains only first order derivatives of 
the respective fundamental solutions (cf. [8] of this paper or the next section 
of this paper). Then derivatives of first order of the fundamental solution 
have Gaussian a priori estimates. We can write down estimates which are 
convolutions, and then we may use the generalized Young inequality in order 



(3 



to establish L 2 estimates for these convolutions (using L 1 estmates of the 
Gaussian). Note that locally we may use the bound 

(26) 



(t - s)H(x - y)n+l-2fJ, 

for 0.5 < n < 1 for the first derivative of the fundamental solution. From the 
classical representations of the functions b~u p,k ' 1 we get a contration property 
with respect to the if 2 -norm uniformly r. This implies that the local limit 
function (fc oo is Holder, and this implies that the limit function has a 
representations in terms of a fundamental solution. Higher order regularity 
may then be obtained by Gaussian estimates and the generalized Young 
inequality using partial integration and the adjoint. □ 

Now in addition we may prove that u P ' 1 are differentiable with respect to 
T. Furthermore we need a bound on the first order coefficients u P,l ~ l which 
is independent of the time step number I. Indeed we have 

Lemma 2.2. 

K' /_1 (/ -l,.)\ a <C for alll<i<n (27) 
for a constant C > which is independent of the time step number I > 1. 

Proof. Again we consider the case n = 3 We observe that the growth of 
w ' 1 is linear with respect to the time step number I. At each time step a 
time step size of order p\ ~ j then ensures that the first order coefficients 
u P,l ~ l {l — 1, .) for the equations for u?' 1 ' have a uniform bound (independent 
of the time step number /). □ 

Invoking classical regularity results for linear parabolic equations and 
induction over / yields 

Theorem 2.3. The Cauchy problem ([TJ) on [0, oo) x W 1 with initial data 
hi € H s for all 1 < i < n and all s € R has a global regular solution 

iii £ C°°([0,oo) x IT) (28) 

along with u(t, .) € H s for all t € [0, oo). 

This is not a new result but the proof is more elementary compared to 
an alternative method which establishes first the estimate 

^\\u(t,.)\\ H s <\\u(t,.)\\ H s + iY, Yl II^ Q ^^IIl2 -2||V U ||^ (29) 

ij \a\+\/3\<s 

At first glance from the construction it seems that the estimate for the 
solution function u increases linearly on a time scale Ylk=i Pk with Pfe ~ r, 
and, hence, quadratically on an uniform time scale. However, at each time 
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step, if we have constructed u p '\t, .) £ H s for 1 < i < n uniformly in r with 
r € [I — 1, 1], then the construction gives 

K'\r,.)\ H u <\uf-\r,.)\ Hk +C k (30) 

with a constant depending on k according to the level of differentiability 
we want to obtain (accordingly the stepsize pi has to be chosen for each given 
k). However, the growth of size at each time step can be compensated 
by a choice p\ ~ jq-. Furthermore, we can easily deduce the existence of 
classical solutions u p ' 1 of (|19p at each time step. First we choose k > m + hn 
in the construction above and obtain spatial differentiability of order m 
uniformly with respect to r. Then the product estimate 

\fg\H-<C a \f\ H ,\g\H. for s> l -n (31) 

gives 

-fr-Pi»2^-^r-pi2^ u i ~^-^ H ■ ( 32 ) 

j=l 3 3 =1 J 

Hence for choice k > 2 + \n we have differentiability of u?' 1 with respect to 
time. The existence of a fundamental solution of the equation 



dpP ,i f sy " ,app.' 



is ensured in terms of the known functions u?' which are indeed now known 
to be Holder continuous with respect to the time argument r and the spatial 
argument x. Hence we have classical representations of the solution and an 
ordinary maximum principle for linear parabolic equations tells us that the 
maximum of the solution function u p ' 1 does not increase over the time inter- 
val [I — 1, 1]. As the step size compensates the linear growth of the estimate 
immediately linked to the original construction this holds independently of 
the time step number /. We note 

Corollary 2.4. (Same assumptions as in theorem \2.3\ above). 
The solution u of f7]) is bounded. 

3 Extension of the scheme to the incompressible 
Navier Stokes equation 

Next we look at the relation to the incompressible Navier-Stokes equation. 
We write it in the same scheme frame as above with time step number I > 1 
and time step size p\ ~ \ at each time step number /. The velocity is 
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denoted by v = {v\, ■ ■ ■ ,v n ) T (recall that T is for 'transposed'). We have 
an additional scalar function p, and the equation system scheme becomes 

dr — Pl v l^j=\ Q X 2 Pi 2^j=l v j dx 3 Pi dXi ' 
diw^ = ( 34 ) 

I v^(/-l,.)=v^" 1 (/-l,.), 

where for each I > 1 we consider this system on [I — 1, Z] x M n . For Z = 1 
the initial conditions are v p,1 (0, .) = v p ''~ 1 (Z — 1,.) = h(.) with the same 
h € H s (M. n ) for all s € M as in the preceding section. The Leray projection 
form of these equations is obtained time-step by time-step by elimination 
of p p ' 1 via an equation for div v p '' which simplifies by incompressibility (at 
each time step I > 1). First, this leads to the Poisson equation for p p,l (r, .) 
for r G [Z — 1, 1], i.e., 

ApP.' = - Y (35) 
^ dxn dxi ' 

where we suppress the notation of the evaluation at r (r serves as a parame- 
ter). Hence at each time step we have the local equation in Leray projection 
form of the incompressible Navier-Stokes equation. If 




dvf dvf 



dxj dxi 



(r, .) (36) 



is in C 1 (R n ) n L 1 (R n ) for all r € [1-1,1], then this is justified by a well- 
known result. Indeed we have 

Proposition 3.1. Assume that n > 3 and that f G C 1 (R n ) n L 1 (R n ) and 
let K be the fundamental solution of the Laplacian, i.e., of the equation 

AK = 5, (37) 

where 5 denotes the Dirac distribution. Then for n > 3 

I ™|2 — n 

K(x) = - ' 1 . — (38) 
v ; (2-n)w„ 1 ; 

(with uj n being the surface area of the unit sphere) defines a distributive 
solution and determines a classical solution w G C 2 defined by 

w(x) := I f(y)K(x - y)dy (39) 
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of the Poisson equation 

Aw = f. (40) 

Moreover, the gradient 

Vw(x) := [ f(y)VK(x - y)dy (41) 

is well-defined. 

Note that for application in our scheme L 2 -theory for the data is enough 
in order to apply the proposition 13.11 since 

for any r € [1 — 1,1]. This means that v p ' E H 1 implies that the left side of 
(|42p is in L 1 . Outside a ball we may then use that the first derivatives of the 
kernel K j are in L 2 in order to estimate pressure terms. The truncation of 
the kernel inside a ball convoluted with the data is easier to handle. Indeed 
since the data are in our scheme at time step I > 1, i.e., v p ' (I — 1, .), are in 
H 2 f~)C 2 , and we have 

" I dv p ' 1 dvf " | fluf 

E|^-V)^-1,)|<^E|^-V)| (43) 

We may then use that the right side of (|43|) is in L 2 and together with the fact 
that the truncated kernel l^ r f \K t i is in L 1 the generalized Young inequality 
leads to the conclusion that the convolution with truncated kernel is also 
in L? . Similar for first derivatives. Next we extend this standard lemma 
in order to meet these purposes. It is useful to have the gradient of the 
pressure in L 2 (resp. in H 1 ) since estimates in terms of convolutions of 
the modulus of the gradient with a Gaussian a priori estimate help us to 
conclude that the convolution bound is itself in L 2 (resp. in H 1 ). Note that 
the relationships in (|14p and (j 15|) make L 2 -estimates for the gradient of the 
pressure useful in order to get L 2 -estimates for the approximations v P ' ' of 
our scheme and their (weak) derivatives up to second order. Furthermore, if 
we want L°°-estimates for approximations v P ' k ' 1 and their (weak) derivatives 
up to second order, then L 2 - estimates for the gradient of the pressure are 
also useful since the Gaussian is in L 2 . We prove 

Lemma 3.2. For a multivariate function f we say that f S H k = H k (W 1 ) 
if all multivariate derivatives D%f with multiindex a = (a±, ■ ■ ■ ,a n ), oti > 
and \a\ = ^27=1 a i — ^ are * n L 2 ■ Similarly, for a multivariate function f 
we say that f € H k ' 1 = H ki1 (W 1 ) if all multivariate derivatives D®f with 
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multiindex a = (a±,--- ,a n ), on > and \a\ = Y17=l a * — ^ are ^ n ■ 
Accordingly we set H° (M. n ) = I? . For k > assume that fa, gj G H k+2 n C 2 
for 1 < i < n, and define for all x G ]R n 



T/ien toe /iaue ^-p / ' 9 G n C 1 and p^ 9 G i? fc+2 ("I C 2 (defined accord- 

ingly). Moreover, 



d f,9 
OXi 



„ k+1 <C + C(\f\ Hk+2 + \g\ Hk+2 ) (45) 



for some constant C > 0. Especially, in case n = 3 and if f,g G i? 2 D C 2 , 
then we have -£-p f ' 9 G H 1 n C 1 and £ H 2 t~)C 2 . 

Proof. We consider the case /, o G 2 H C 2 and n = 3 which is essential. 
Note that the second derivatives of / and g are in L 2 n C, and therefore in 
the closed space of continuous functions vanishing at infinity. Hence they 
are bounded. The basic idea is is the following. The representation of the 
gradient of the pressure is a convolution of a (sum of) products of 

E lr¥ < 46 > 

^— ' OXj OXi 

1,3=1 J 

with the gradient of the kernel K. Since for all x G R n we have 



e- j=1 Stfe(-)<E- fc=1 HSj"w + Hiy"(x), (47) 

it follows that the (sum of) product(s) function in (|46p is in H > and its 
first order derivatives are in L 1 (since fj and gu are in H 2 ). In order to have 
the gradient of the pressure in H 1 the next idea is to split up the kernel K 
(or a first order derivative of it) into two summands writing the gradient of 
the pressure functional as 

ity- 1 ^ = >::;.<■ . (Ifcft) o* - v^m* 

= k« E], k =i (StSj) (* ~ y)My)KAy)dy (48) 



2 



along with a smooth function supported on a ball of radius r > around 
B r (0) and zero elsewhere. Recall the Einstein notation where K i denotes 
the partial ith derivative of the kernel function K. Note that the truncated 
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kernel 4>B r K^ is in L 1 but not in L 2 . Concerning the L 2 -estimate of the 
gradient of the pressure functional this is no problem. For dimension n = 3 
we have that (j>B r K is in L 2 . Therefore, using convolution rule we may write 
the first term of (1481) in the form 



ELl ((Ifeft) {*-y)M{y)) K(y)dy. (49) 



We may then use local L 2 -integrability of K = K% (dimension n = 3), since 
for the truncated integral we have 

/ K 3 (x)dx~ [ \s\ i - 2n s 2 ds~r. (50) 

J B r JO 

This is of use if we expand YIj=i ELi ((J^rffj) ( x ~ y)^B r ){y)) . and 
write (|49p as a sum of convolutions defining an appropriate function (f>B r 
which has a compact support of first derivatives. If we consider derivatives 
of the gradient of the pressure functions we need another idea, but this works 
for the estimates within the gradient itself. Well, the localized terms are not 
really a problem since we can invoke local regularity theory for the Poisson 
equation. However, since second derivatives of the source function fj,gk are 
continuous, there is also another possibility which also works for the second 
derivatives of the pressure functional p^ 9 . We may write 

3£ = kn EU ELi ((St§f ) - f)) M(y)K,(y)dy 

(51) 

< I^nZ n k=1 c\^(x -y)\\4> Br )(y)K t (y)\dy, 

for some C > 0, and where we use the boundedness of / £ H 2 n C 2 and 
its derivatives up to second order. Now the left side of (j63j) is in L 2 by 
the generalized Young inequality, since Ki G L 1 and second order partial 
derivatives of gu are in L 2 . The second term on the right side of (|48jl may be 
estimated using the fact that the function (1 — <f)B r )K,i is square integrable 
since for n = 3 we have 

i 1 - M-r-k^dx ~ -FS 2 ds <oo. (52) 

n \B r (0) Fl Js>r s 

Now let us fill in some more details of this argument. First we specify 4>B r 
to be 4> t where we write 

= k« Ei,fc=i (ftft) - v)&W*(y)dv (53) 
+ 4* E? lfc =i (gfefe) (x - y)(i - Mv))KM*v, 



12 



where 



Note that we have 



f 1 if \y\ < e, 
else. 



if e 2 < \y\ 2 < 2e 2 



(54) 



(55) 



because , gu € H 2 n C 2 and 



En 
7,fc=l 



(x) 



(56) 



From proposition 13, II we conclude that p-^' 9 € C 2 , hence -§^r.V^ 9 € C 1 . There 
are different methods to estimate the first summand on the right side of (|53|) . 
We have seen one above. Another variation of argument is the following. 
We write 



IEIh (Hit) (x-y)Uy)KAy)dy 

= / R » Elk=i (ffiSj) 1 {i/ Jffc c*-y)i>i & ijyNiy}* 1 " y)<t>e{y)K,i{ y )dy 

+ /r« E^=i (Stlfj) 1 {i/ J , fe (x-j / )i<i or ia w (as- 1 ,)i<i}( a; ~ y)My) K Ay) d y 

(57) 

?or the first term on the right side of ()57|) we get for small e > 

I Ei, fc =i (ftfe) ^i/^-^i & i^-^i}^ - y)Uy)K,i(y)dybi 9 \ 



(58) 

for some compact set K C W 1 because fj , € if 2 D C 2 . Since fj , g% G C 2 
we have 



dfj_dgk\ 1 ^2 
dxkdxjj 

With the same argument we get that the latter function in ii 1 (differentiate 
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the product and proceed as before). For the second term of (|57p we have 



Jr« Ylj,k=l 
— Jr™ Ej',fc=l 



dxk dxj 



H\f,^-y)\^ or | aw (x-»)|<i}( a? -l')^(f) 

{ii 9feJ (^)i<i}( x - y)^(y)|^(y)|^ 



j,fc=l 



1 {l/M(*-»)l<l} (iC ~ I/ ^ 6(!/) 



dy 

dy. 



(60) 

Hence since 4> t K^ G L , and |<7fej| G ^ 2 an d G L 2 the function corre- 
sponding to the term in on the left side of (|60p is also in L 2 by the generalized 
Young inequality. Again an analogous argument shows that this can be done 
for the second derivatives in the analysis of the pressure functional p p <f' 9 . 
As we said, there are variations of this argument. Next to the arguments 
given for dimension n = 3 we may use 



dfL§2k 
dxk dxj 



En 
j,k=l 

< cc 2 Elk=i 



L 2 



< C Elk=i fj9k 



11, 



(61) 



by a standard estimate 



ffj 



H" 



: C 



9k 



h 2 



for s > \n This may be 



used for estimation of the first term on the right side of [53] were we have 
4> e K i G L . Hence, we have 



x — > p 



E 

j,k=i 



dfj dg k 
dxk dxj 



{x-y)4) e {y)K^(y)dy G L 2 



(62) 



by the generalized Young inequality. However, for second derivatives we may 
need to describe the scheme in H m for larger m in order to get estimates for 
the second derivatives. Note that we may also apply local regularity theory 
for Poisson equations for this part. 

For the second summand term of the right side of (|53p we may write 
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with our specific function <p e 
Jr3 E"=i ELi {fj,k9kSi, m 0* ~ ~ ( l ) e)(y)K ! i(y)dy 

Jr3 E"=i ELi (fj,k,m9k,j + fj,k9k,j,m) - y)(l - 4>e)(y)K ti (y)dy 

< E"=i ELi I (& lB1 + sgj + /J* + sk, m ) - »)|(i - *«)(v)#,*(v) 

(63) 

Then we may use the generalized Young inequality (this time observing 
that |(l-0 e )(y)^| € T 2 and ££ =1 ELi (/ilj) (•) e ^ and we conclude 
that the function corresponding to the left side of (|63|) is in L 2 . Again 
their are variation of this argument where me invoke second derivatives of 
(1 — (p t )(y)K t i. Summing up we get an upper bound of form 

C + C[ max \gi\ H 2 + \fj\ H2 j (64) 

for some generic constant C > as desired. The proof for fcgj € H k+2 with 
k > is also similar. □ 

We need also standard Gaussian estimate for the fundamental solution 
and its first derivatives. We have 

Lemma 3.3. Let D := [To, T\] x O C R n be a domain along with T\ > To > 
0, and let 

_ dp sr^ d z p s-^ dp 

^ = ^-£^-^-+E^ = ° ( 65 ) 

be an equation which satisfies 

i) L is uniformly parabolic on the whole of [To,Ti] x Vt, 

ii ) the coefficient functions aij are uniformly Holder continuous with Holder 
constant a/2 € (0,0.5) with respect to time and Holder constant a G 
(0, 1) with respect to the spatial variables, i.e., a^ & (ja/2,a 

Hi) the coefficient functions are Holder continuous with Holder constant 
a € (0, 1) and uniformly with respect to time t. 

Then a fundamental solution p of equation 165\) exists and satisfies the Gaus- 
sian a priori estimates 

C ( X(x-y) 2 \ , . 

\p{t,x;s,y)\ < — — Tjexp — — — , (66) 



dy. 
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and 



_d_ 



p(t,x;s,y) 



< 



C 



(t - s)(™ +1 )/ 2 



exp 



4(t - s) 



(67) 



for some constants C > and A > (X less or equal to the lower ellipticity 
constant in general). Note that for t > s these a priori bounds (as functions 
of x — y) with have a uniform bound in L 1 . 

We also use some related standard results on the adjoint. We have 

Lemma 3.4. Assume that the assumptions of the preceding lemma are sat- 
isfied on a domain D = [Tq,Ti] x M n ; and that in addition we have 

i) The coefficient functions aij and their first and second derivatives are 
bounded continuous functions on D, 

ii) the coefficient functions b{ and their first and second derivatives are 
bounded continuous functions on D. 

Then the fundamental solution p* of the adjoint equation L*p* = exists. 
Furthermore, for t > s 

p*(s, y; t, x) = p(t, x; s, y) (68) 



and analogous relations for the partial derivatives hold. Moreover, and sat- 
isfies the Gaussian a priori estimates 



\p*(t,x;s,y)\ < 



C 



(t 



and 



-^-p*(t,x;s,y) 



< 



W2 



c 



exp 



(t - s)^ 1 )/ 2 



exp 



A(x - y)' 
4(t - s) 

4(t - s) 



(69) 



(70) 



for some constants C > and A > (X less or equal to the lower ellipticity 
constant in general). 

Concerning these lemmas we note that in our scheme only the first order 
coefficients are variable. The useful relationship of the fundamental solution 
and its adjoint may also be verified directly using the Levy expansion as we 
did in [8] . We also discussed the matter in [7] . Our lemma 13.21 showed that 
we may base a global scheme for the incompressible Navier-Stokes equation 
on approximations in Leray projection form which are H 2 with respect to 
the spatial variables. The scheme is in H 2 for fixed time with respect to 
the spatial variables (which is easier to see), and this helps us in order to 
estimate the gradient of the pressure functions of type H as in the 

lemma 13.21 above. 
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Next, for each I > 1 we consider the Lery projection form of the incom- 
pressible Navier-Stokes equation at each time step I > 1, i.e., 



( dv?' \-^n d 2 v p i ' 1 ^ n 

-tfr = pi v 2^3=1 ~d£- ~ Pi 2-j=i 



p,l dv 



j dxj 



^P 1 JlR n Ei,i=l 



(71) 



k v^(Z-l,.) = v" ,_1 (Z -!,.)■ 



Next we may define a scheme similar to the scheme of the multivariate 
Burgers equation and with time step size p\ ~ j . Having computed v?' 1 — 
1, .) G H 2 fl C 2 we set up a local iteration scheme. In general, for n > 3 a 
condition of form i^'' -1 (Z — 1, .) G fP 1 n C 2 for an integer m with m> \n 
is an appropriate choice. But this is quite similar. Hence, think of n = 3 
and of data and functional series approximations with u?' ' (t.) € H 2 first. 
Again, the solution in Leray-projection form is constructed via a functional 
series 



p.i p 
v- = v 



n, 



fc=i 



where uf' 1 ' solves 



(72) 



p,l,i 



2„,P,1,1 



9 2 « 



^ Jr™ Ei,j=l 



(l-l )y )£-K n (x-y)dy, 



k W-'(Z-l I .)=v , - 1 (i-l,.) ) 



and owp ' = uj' — ' , 1 < i < n solves 



d&v 



p,k+i,i 



Pi [ 2.7=1 a^f 



En p,k,l d&v. 



p,k+l,l 



f).r: , 



-piE 7 - 



,p,k,l dvP*- 1 ' 1 
dxi 



Pi J R n i^n,i(a; - y) ( (Ej, fc =i + <'j M ) ( r - y)) x 



[ j v p.*+i.«(Z_ i,.) 



(73) 



(74) 



17 



and where Svf ,k+1 ' = vf' k+1,1 — v^ k)l for k > 1, and 5v P ' 1,1 = v?' 1 ' — v p ' 0,1 := 
1 it j j 

Vj' 1 ' — v P ' ~ (I — 1, .). Note that Sv^' 1, = v?' ' X — hj at the first time-step. 

The equations for v?' 1 ' 1 and for 5v p,k ' 1 are linearized and localized equations 
where by localisation we mean the fact that the global integral terms in 
the equation for v?' 1 ' and for 5v?' ' are given in terms of the initial data 
and the data from the previous iteration step respectively. Assuming that 
v i' l ~ 1 (l ~ 1) •) £ H 2 we show that v p,1 ' 1 (t,.) in H 2 uniformly with respect 
to r and that the series (^5v p ' k ' l (T, .)J satisfies a contraction property in H 2 

uniformly with respect to r € [I — 1,1]. We show that the functions v?' 1,1 
are limits of a functional series where we have a contraction property for 
the elements of the series of form Sv p ' k+1,1 . This implies that the series (|72p 
evaluated a r € [I — 1, 1] converges in H 2 such that for n = 3 we may apply 
the following extension of a standard Sobolev lemma. 

Lemma 3.5. For s = a + k + \n with a £ (0, 1) we have 

H s c C a , (75) 
where C a is the space of Holder- continuous functions. 

We note that the IP estimates which are useful for higher dimensions 
may be used in the context of a more general lemma. We have 

Lemma 3.6. For s > k + ^ with a G (0, 1) we have 

H s > p C C k , (76) 

where H s,v is the space of functions f where A s f S LP along with 

A s = [/-2vr- 2 A] s/2 . (77) 

Applying such type of lemmas we can ensure that the first order coeffi- 
cients (evaluated at time r) of the equations which determine our approxi- 
mations v? ,k ' 1 satisfy classical conditions which are sufficient for the existence 
of fundamental solutions of the associated linear parabolic equations of our 
scheme. Especially first order coefficients evaluated at r are in the Holder 
space C a for some a € (0, 0.5) and in Qj uniformly with respect to r, i.e., the 
in the space of continuously differentiable functions which vanish at spatial 
infinity. 

The construction here defines a weak solution in H 2 and in H 2 '°° as 
n = 3. We have not mentioned ff 2,00 -estimates explicitly, but the estimated 
above can be adapted straightforwardly. Note that the application of the 
generalized Young inequality is even more simple in this case: for r = oo we 
have l + f = | + |, and we may use p = q = 2 outside a ball and p = oo 
and q = 1 inside a ball. Since v p ' 1 £ C a as the limit of the functional series 



IS 



in (|72p , we have representations in terms of v p ' 1 in terms of the fundamental 
solution of 



dp I A d' 2 pi -A i dp 

and this leads to the immediate conclusion that the solution is classical. 
Note, however, that for local restrictions to a bounded domain f2 C R"' the 
series v p ' k,l (r, converges to a limit v p ' 1 (t, in a classical Banach space. 
Recall the following fact. 

Proposition 3.7. For open and bounded Q C R" and consider the function 
space 

C m (O) := |/ : ft -> R| <9 a / exists for \a\ < m 

(79) 

and d a f has an continuous extension to $1 j 

where a = (ai, • • • , a n ) denotes a multiindex and d a denote partial deriva- 
tives with respect to this multiindex. Then the function space C m (JTj with 
the norm 

|/U := l/lo»(n) : = E l^^l ( 8 °) 

a|<m. 

is a Banach space. Here, 

\f\ :=sup|/(x)|. (81) 

This leads to a second argument that the limit is indeed of form v?' (r, .) G 
H 2 nC 2 uniformly in r and such that v p ' 1 satisfies the incompressible Navier- 
Stokes equation locally on [I — 1,1] x R n . This variation of argument has 
the advantage that it does not depend on dimension. On the other hand we 
can do the construction in H m instead of H 2 for m > 2 + ^n, so this is a 
matter of taste. The last step then is to show that we have a linear bound 
of growth with respect to H 2 . We shall show that 

\v p /(l,.)\ H , <\v p /-\l-l,.)\ H2 + C 2 (82) 

for a constant C2 which is independent of the time step number I and which 
holds for all 1 < i < n. This implies that we have a global bound 

max \v P ' l (l, .)\ H 2 < max \hi\ H 2 + IC2 (83) 

l<j<n l<i<n 

Now let us consider this program of proof in more detail. First we observe 

Lemma 3.8. Let vf~ l (l - 1, .) € H 2 n C 2 for all 1 < i < n. Then there 
exists a classical solution v p ' 1 of FHfy with 

/• y (r,.)eff 2 nC 2 (84) 
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for all t £ [1 — 1,1]. Moreover, for n = 3 we have 

v p,1,1 {t, .) G H 2 '°°, (85) 

where the latter space denotes the Sobolev space with weak derivatives up to 
second order in L°° . 

Proof. Since v r i ' p ' l ~ 1 (l - 1, .) G H 2 we have < W_1 (£ - 1, .) G H 2 n C Q for 
a £ (0, 0.5). Hence the fundamental solution p l of 

j=l 3 j=l J 

exists (constructible in the classical sense by the Levy expansion) and the 
solution of the Cauchy problem in (|73p has the representation 

v? x \t,x) = / R «wf i_1 (Z - l,y)p l {r,x,l - l,y)dy 

(fi„pj— 1 a p,i— i \ 
I fe % j(f-i,y)gn,i( g -y)x (87) 

xp l (s, x, I — 1, z)dzdyds. 

Here, recall that iT ni j denotes the partial first order derivative of the kernel 
K n with respect to the ith variable. Hence we have uf" (r,.) G C 2 for all 
r G [/ — 1,Z] which follows from classical analysis of the Levy expansion of 
the fundamental solution where we may differentiate under the integral in 
order to get a representation for the derivatives of first order for r > I — 1. 
Moreover, the second derivatives of the last integral in (|87|) have an adjoint 
representation (cf. also the argument in [?] and [7]) such that the second 
derivatives of u?' 1 ' with respect to the spatial variables and x m is the 
sum of 

/ ^'^-^f^T^x.l-l^dy (88) 
(which exists since vf' i ~ 1 (^ — 1, ■) is Holder), and the second summand 

+ Pi Ii-i Jr« Imp- 52p,j=i TEl i^T^j ^ ~~ ^v)~5^ K n{ z ~y) x 

x Q§—p l '*(s, x, I — 1, z)dzdyds 

(89) 

(for the adjoint p''* cf. also [8] and [7] ). The term 
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corresponds to an L 2 -function according to our lemma 13.21 above (and to 
a L 2 -function if n = 3), and Gaussian estimates for the first derivatives 
of the fundamental solution and its adjoint plus an application of Young's 
inequality in its generalized form ensure that (|89p is in 1? for each r (first 
the integrand is in L 2 for each s £ [I — 1, r] and then the integral up to r 
is in 1? where r is considered as an parameter). Note that the Gaussian a 
priori estimate of the fundamental solution in lemma is L 1 C\L 2 for fixed 
r > s. Similar for the Gaussian a priori estimate in lemma l3T4"l Let us look 
at the second term more precisely since this is the term which defines the 
extension of our scheme for the multivariate Burgers equation. We have 

+ pi II-i Jk« Jr« mi Sij=i (~ fe~ J (l ~ ^y)~s^ K n( z ~y) x 



dx 



? p l '*(s,x,l - l,z)dzdyds 



x (t _ B)( n + i)/2 exp (~^ff ) |^(is. 

(91) 

Now our lemma 13.21 and the generalized Young inequality in (|14p and (|15p 
with r = 2 and p = 2 and g = 1, i.e., 

/e^andffGL 1 -> / * 5 G L 2 , if 1 + i = 1 + -, (92) 

and 

\f*9\v<\f\v>\g\v- (93) 

where / corresponds to the Leray projection solution of the gradient of the 
pressure analyzed in lemma 13.21 and g corresponds to the Gaussian a priori 
bound where we observe that for t > s we have 

for p > 1. Similarly for the first term. Moreover, from the representation 
of both summands we observe that v?' 1 ' G C 2 . We note that for n = 3 we 
may apply lemma 13.21 and the generalized Young inequality in (|14|) and (|15|) 
for r = oo and p = q = \ in order to get L°° estimates. □ 



Next in order to construct a local solution of the Navier-Stokes equation 
we establish a contaction property for the correction functionals 5v1' k ' 1 of 
the first linear approximation v^ 1 ' 1 considered above. 
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Lemma 3.9. If 1 (l - 1, .) G H 2 n C 2 for all 1 < i < n then there is a 
contraction property 

ISv^^.^K^Sv^iT,.)]^. (95) 

Moreover, if ^f' (Z - 1,.) G # 2 n C 2 Zfcia is a contraction in |.|2 suc/t 
t/iat restrictions of the functional series to a bounded domain f2 converge 
in a classical Banach space of twice differ entiable functions with continuous 
extension at the boundary , i.e. ^uf' fc, '|n^ converges in the classical Banach 
space described in proposition \3. 7| for m = 2. Moreover, for n = 3 we also 
have the contraction estimate 

\5v^\r,.)\ H ^ < 1 -\5v^ 1 {t,.)\ h ^. (96) 

Proof. We prove the theorem in case n = 3. This implies that we can apply 
the product rule f)31[) with s = 2 > | = ^. The generalisation to dimension 
n > 3 can be done similarly but with fP ,p -estimates (cf. our remarks on 
.fP' p -estimates and lemma [3721 above) . At each substep k of time-step I > 1 
we may assume inductively that 



vf k \ r , .) = vf' l (r, .) + $v p > m > l (T, .) G H 2 C C a (97) 



m=2 



for a G (0,0.5) and uniformly with respect to r G [/ — 1,1]. Moreover, we 
know inductively that Vj ,k ' l (T, .) G C 2 for all r G [Z — 1, Z] . For k = 1 (when 
we interpret the second summand in (|97|) to be zero) we know this from the 
previous lemma. Hence inductively with respect to the subiteration index k 
we know that the fundamental solution p k ' 1 of 

j=l j=l j 

exists, and it follows that the solution of the linear problem (|74p has the 
representation 



x 



-pi JJLi J" R » Ej (^i' fc ' z ^fe— J ( s > y> fe, '(^ x, s > 

(99) 

Jj-i /r» J r « - y) ( (Ei, fc =i (<j ' + v kTj ( T > y ">) x 



( Sv j'k' l ( T > yj) )p k ' l ( T > x > s > z)dydsdydz. 
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For the first order derivatives we have the representation 

-pi /,li S R n Ej ( Sv j ,k ' l dv 'a Xj ' j ( s , y)di^p k ' l ( T , x , s , y) d y ds 

+PI /,ll / H » / R » #n,i(* - y) ( (Elp=l ( V k k p + ( T > V)) X 

(6vf *' l (T, yfj ) Q§^P k,l (T, x, s, z)dydsdydz, 



(100) 



and for the second order derivatives we have the representation 
^5vf +1 \r,x) = 

+Pl JT-i /«« Ej ot {*vf' lQ ^) v)^^' S , y )dyci S 

/ill /r« /k» ^n,m(^ " (E£ P =1 (<f + ( T ' f) 

(6v?'p' l (T, yj\ ) -£^P k ' l '*(T, x, s, z)dydsdydz, 

(101) 

and where p k ' 1 '* denotes the adjoint (consider also part I of this article). The 
following contraction estimate involves an estimate for 



p,k,l p,',l 

vy = vy 



+ Yl 6v i ,m,L ( 102 ) 



m=l 



because of the nonlinearity in the first order coefficients. This involves the 
estimation of v\ ,p ' 1 . Recall that we have the representation in (|87p such that 



v-' 1,1 (t,x) - v p '°' 1 (t,x) = f Rn v?' 1 1 (l - l,y)p l {r,x,l - l,y)dy - v p '°' L (t,x) 

-fci — krj ( l - !. y)ax- K n(z -y)x 

xp l (s, x, I — 1, z)dzdyds. 

(103) 

Well, we defined v p '°> 1 (t, x) = v p,l ~ l (l — 1,.). We could have defined v?' 0, 
such that the first summand on the right side of ()103p cancels. Anyway 
classical analysis tells us hat there is a bound in the relevant norms. The 
essential term is the second summand which we may estimate using lemma 
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13.21 For the equation in (|104p we get the estimate 



d 2 r P ,k+l,l, V 

x m dx q 0V i 

^ Pi Ii-i Ej a 



L 2 



9_ / ^ v P^ l ^± 



p,k-l,l 



<).l- j 



■(«,-) 



( r _ s )(n+l)/2 ex P ^ 4(r-s) 



(is 



( T _ s )(n+i)/2 ex P ^ 4(r-s) 



p,fe-l,i 



L 2 



CoC'ds 



+pi IU | / K n - y) (^- (EI p =i (<f + <J M ) (r, y)) x 



<Pz 



L 2 







p,k-l,l 
v i 




Jf 2 





CoC 



p,k,l , p,k—l,l 

Vp ' +vp 



5v p 



k.l 



H 2 



CC 2 C 



Vp ' +v!p' 



1,1 




$ v P,k,l 


c* 




H 2 




H 2 



<p;E r 

(104) 

where we use the product rule (|3ip above with the constant C2, the uniform 
bound 

C ( A(.) 2 



( r _ s )(n+i)/2 exp ^ 4(r-s) 

lemma [3721 (therefore an additional constant C for the second summand), 
and the generalized Young inequality in (|14|) and (|15|) above. Moreover, we 
may choose C* > CC2C 1 '. For generic C* if follows that 



(105) 



p,fc+l,Z , 



ff2 



W 2 



k.l 



H 2 



Note also that we used that 



-kA t1 <c 



(106) 



(107) 



for some constant C > (this does not hold for L 2 , but leads to L 2 -estimate 
for a convolution with convoluted L 2 -data via the generalized Young in- 
equality). Next, in order to estimate PiC*^2 



p.k.l 1 p,k—l,l 



we first 
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estimate 



p,l,l . p,0,l 



H 2 



H 2 



E P v£ 1 > l -v p p > l - 1 (l-l,.)+2v p /~ 1 (l-l,.) 



(108) 



and then proceed by induction over k showing that a certain choice of pi is 
appropriate. We assume inductively that for all 1 < i < n 



p,Z — 1 



(l-l,.)\ H 2 <C* + (Z-1)C* 



(109) 



for some generic constant C* > and Z > 1. Note that for I = 1 we have 
|/ij(Z — 1, -)\h2 < C* for all 1 < 2 < n for some generic C*. In order to start 
the (sub-)induction for k = 1 we go back to representation in (|103|) . We 
have to estimate 



< 



H 2 



1 (l - l,y)p l (r,x,l - l,y)dy - v p ' 1 1 {l 



dxj dxi 



(l-l, y )d-K n (z-y)x 



H 2 



vf\r,.)-v^ l {r,.) 
+ Pi Il-l Jr« /l" 12i,j=l 

xp l (s, .,1 — 1, z)dzdyds 

(110) 

for all r G [1 — 1,1]. For the second term we may use the techniques of 
lemma l3^2t especially the .ff 2 -estimates of convolutions together with a priori 
estimates of the fundamental solutions of adjoints and their first derivatives 
in L 1 -norm. This together with the generalized Young inequality leads to the 
estimate of linear growth. Note that we involve i? 2 -estimates since we have 
linear growth of the function v p,l ~ l (l — 1, .) with respect to the H 2 -norm. 
Especially, second derivatives are represented and estimated such that the 
highest order terms (partial derivatives of second order with respect to the 
variables x m and x q of the |.|^2-norm become 



Pi Il-l J*K n Im n Sj,p=l dxi ( dx„ cL, ) 1>2/) d Xi K n,m{z y) x 



xp\q{s, .,1 — 1, z)dzdyds 
Pi Ii-i /r« Jr™ Ylj, P =i 



L 2 



< 







2 

+ 






dxpdxi 





(*-«) 



(n + l)/2 eX P 



■^^)dzdyd. 



L- 



( l - 1 ,y)-B x - K n,m( z -y) x 



(111) 
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The latter term is a double convolution where we may estimate the Gaussian 
a priori estimate in an L 2 -norm and the convolution function function 

-l „ 




+ 



dv p v 



1-1 



dxj 



(!-l,y)-£rK ntm (z-y)dy (112) 



in L 2 -norm along the techniques of lemma 13.21 invoking the generalized 
Young inequality in order to get an L 1 -bound n 2 C(C* + (I — 1)C*). Note 
that this is not trivial because K n>m (.) is only locally in L 2 , i.e., is in L 2 QC , 



so this requires the techniques of lemma 13.21 Next note that the first or- 
der coefficient of the equation for the fundamental solution p l are of form 
/9;f['' _1 (Z — 1, .) which has a i7 2 -bound independent of the time-step number 
I for the choice pi ~ j . Then applying the generalized Young inequality once 
again we get (with another generic C-factor from L 2 -integral of the Gaussian 
and collecting summands of the H 2 -norm) the inequality 

Pt Ii-i Im n /r« Ylj >p =i 



dx v 



du 



P ,l-i 



(l-l, y )^-K n (z-y)x 



(113) 



xp l '*(s, .,1-1, z)dzdyds < n 2 C 2 (C* + (I - 1)C 



H 2 



f Rn v?> l -\l - l,y)p l (r,x,l- l,y)dy - vP' l ~\l - 1, .) 



Furthermore, looking at the first summand on the right side of (|114|) we find 

< C* (114) 

n - 

with generic C* > and choice 

Pl = A(n 3 C 2 C*(C* + IC*)) ^ 15 ^ 

for some constant C* and greater than maxi<j< n \ hi\n2 ( which otherwise 
depends on dimension n and the parameters of the Gaussian a priori esti- 
mate). With this choice we get 

1 



5vf\r,.) 



< 

H 2 ~ 2 



5v p 



i.i 



H 2 



and 



Sv 



from the definition of v p,k ' 1 immediately by induction on k where analogous 
Gaussian estimates involving the fundamental solution p k > 1 and its adjoint 



P,k+i,l 



(r,0 



1 

< - 

H 2 ~ 2 



H 2 



(116) 



(117) 



can be invoked. 



□ 



We note that we may extend this lemma by analogous methods to 
Sobolev spaces which is useful especially for n > 3. The extension is es- 
pecially easy if we consider the H m estimates only at integer time points 
I > 1 which is sufficient if we add a little classical regularity theory for the 
appoximating functionals. 
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Lemma 3.10. Let v?' (I — 1, .) € -fP 1 For I > 1 and pi as chosen above 
we have 

vf\l, .) € H rn and 5v p ' k ' l (l, .) € # m (118) 

and 

\8vf>\l,.)\ Hm <\, (119) 

and for k > 1 

|fef M (/,.)l^<^|^ fc - U a,.)l^ (120) 
/or all 1 < i < n and all k > and m € N. 

We have indeed linear growth on a time scale pi ~ j. 

Lemma 3.11. For all m > there is a a constant C m independent of the 
time step number I such that 

K\l, .)\ H ™ < \vf~\l - 1, .)\ H m + C m (121) 

For m = 2 the proof of this lemma is included in the proof of lemma 
13.91 above. From the latter lemma we conclude that we have linear growth 
on a time scale defined by the time step size p\ ~ j. This implies that 
the first order coefficients of the linearized equations for v? ,0,l ,l < i < n 
are uniformly bounded, i.e., we have for any Holder norm |.| Q with Holder 
coefficient a G (0, 1) 

\p lV P> l -\l-l,.)\ a <C (122) 

for some constant C > independently of the time step number I. Moreover, 
we have Holder continuity of all correction terms 5v? ,k ' 1 , 1 < i < n, k > 1 
uniformly in r € [1 — 1,1] and independently of the time step number I > 1. 
This implies that the scheme is global. 

The argument we have proposed here is considerable simpler then the 
argument in [Uj where a dynamic control function is used. On the other 
hand, the dynamic control function may stabilize the scheme and allows for 
a uniform step size. It may also be used in other situations. Adding a control 
function as in [6] also leads to a globally bounded solution immediately, while 
a little additional argument is needed in the line of argument of this paper. 
Note that the present argument also leads to a different proof of the classical 
proofs for global L 2 -existence in [3] and [9]. It would be interesting to apply 
the scheme using probabilistic methods considered in pQ, [2], [5], and [3]. 

1 kampen@wias-berlin . de , kampenOmathalgorithm . de. 
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